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Recent high precision experimental data for heavy-ion fusion cross sections at en- 
ergies in the vicinity of the Coulomb barrier systematically show that a strikingly 
large surface diffuseness parameter for a Woods-Saxon potential is required in or- 
der to fit the data. We discuss possible origins for this anomaly, including the 
effects of dissipation and the sensitivity of fusion cross sections to the choice of 
inter-nuclear potential. Our study suggests that the frozen density approximation, 
which is often used in analyses of heavy-ion reactions, may have to be re-examined 
for heavy-ion fusion. 



1 Introduction 

One of the challenging problems in heavy-ion physics is to reproduce simul- 
taneously experimental data for several processes, such as elastic and inelas- 
tic scattering, fusion, and Daxticle transfers, within a unified framework. A 
coupled-channels approach Elm provides an ideal tool for this purpose. The 
most important ingredient in the coupled-channels calculations is the inter- 
nuclear potential. A very popular parametrization for this is the Woods-Saxon 
form (WS), given by 

^Ar(r) = -K)/[l + cxp((r-i?o)/a)]. (1) 

This potential has been used widely in analyses of heavy-ion reactions and has 
enjoyed success. 

For scattering processes, it hasJ|ifien well accepted that the surface diffuse- 
ness parameter a is around 0.63 fm.clij In marked contrast, recent higlxprecision 
experimental data for fusion cross sections, which were measujedBlj'u with the 
aim of extracting experimental fusion barrier distributions □, suggest that a 
much larger diffuseness parameter, ranging between 0.8 and 1.4 fm, is needed 
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is is not only for a particular system, but seems to 



in order to fit the 
be a general trend. 

We illustrate this problem iruFig. 1 by comparing experimental data for 
the ^^O -|- 208 fusion reaction eI with various calculations with the WS po- 
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Figure 1: Comparison of potential model with the experimental data for the + ^"^Pb 
fusion reaction. The upper panel is in the linear scale, while the lower panel in the logarithmic 
scale. The solid and the dashed lines are obtained by setting the surface diffuseness parameter 
of Woods-Saxon potential to be a=1.12 and a=0.63 fm, respectively. The experimental data 
are taken from Ref. 8. 



tential model, which inchidcs the incoming wave boundary condition (see Sec. 
2) but no channel couplings. The problem is twofold: at energies well above 
the barrier, where the fusion cross sections a^us is relatively insensitive to the 
couplings, a WS potential with a diffuseness a =0.63 fm significantly overesti- 
mates fusion cross sections (dashed line) . Changing the Vb and i?o parameters 
is not helpful, since it merely leads to an energy shift in the calculated fusion 
cross sections without significantly changing the energy dependence. On the 
other hand, a potential with a— 1.12 fm (full line) fits the data well. Below the 
barrier, the slope of the exponentially falling cross section is given well by the 
o=1.12 fm calculation but is not steep enough in that of a=0.63 fm. A better 
fit to this region would be given by including channel couplings but only if the 
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slope were correct. Thus both above and below the barrier the much larger 
diffuseness parameter a=1.12 fm is required to match experimental data. 

The discrepancy may have a number of causes that include i) channel cou- 
pling effects which were not included in the calculations, ii) the deviation of 
the nuclear potential from the Woods-Saxon parametrization, and iii) friction 
effects at high angular momenta. Of course, a potential with such a large dif- 
fuseness parameter should be regarded as an effective one, which is appropriate 
only for describing the particular channel. But, it is important to understand 
the reasons for the large discrepancy in the diffuseness parameters extracted 
from scattering and from fusion, in order to improve our understanding of 
the fusion process and to allow reliable predictions of subbarricr fusion cross 
section. 

Among the thxee possibilities above, the first point was discussed by Es- 
bensen and Back.Ej They showed that the inclusion of collective excitations 
in both the projectile and the target for the ^^O + ^"''^Sm reaction in the 
coupled-channels calculations significantly improves the fit, leading to a sur- 
face diffuseness of a=0.64 fm. However, the effects of channel coupling on the 
surface diffuseness parameter of the inter-nuclear potential were not shown to 
be significant for other systems, e.g. ^^O -I- ^'^^Sm reaction, as was discussed 
in Ref. 14. 

In this contribution, we study the other two possible explanations of the 
surface diffuseness problem. In order to assess the sensitivity of fusion cross 
sections to the parametrization of the inter-nuclear potential, and to discuss 
the effects of deviations from the Woods-Saxon shape, we_sjdl]_calculate in Sec. 
2 fusion cross [Sections using a double-folding potential taoEJ, as well as an 
error function. til'tj In Sec. 3, we will discuss the effects of friction on heavy-ion 
fusion reactions using the critiral distance model of Glas and Mosellij as well 
as the surface friction model. c3 The summary is given in Sec. 4. 

2 Choice of inter-nuclear potential 

In this section, we use the barrier penetration model to study the sensitivity 
of fusion cross sections to the choice of inter-nuclear potential. Instead of 
introducing an imaginary part to the inter-nuclear potential, in the calculatiDpa 
shown below, we use the so-called incoming wave boundary condition. ac3 
With this boundary condition, it is assumed that there is only an incoming 
wave inside the Coulomb barrier due to a strong absorption, which should be 
reasonable for heavy-ion reactions. In this model, fusion cross sections are 
identified with the absorption cross section. 
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Figure 2: Comparison between the double-folding potential and Woods-Saxon potentials for 
the ^^O + 20Spb reaction. The upper panel shows the total potential for s-wave scattering, 
while the lower panel shows only the nuclear potential. 



2.1 Double- folding potential 

A double-folding potentiallll'El'0 is constructed by folding an effective nucleon- 
nucleon interaction with the target and the projectile densities: 

VNir)^ J dridr2v{r2- ri~ r)pi{ri)p2{r2) . (2) 
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This potential is intimately related to the resonating group method (RGM) 
and has been widely used in microscopic studies of heavy-ion reactions. 

Figure 2 compares the double-folding potential for the ^^O -I- ^"^Pb re- 
action with Woods-Saxon potentials for two different surface diffuseness pa- 
rameters. For the nucleon-nucleon interaction v, we use the clensity-dependent 
Michigan three-range Yukawa (DDM3Y1-Paris) interactions, together with 
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the zero-range approximation for the exchange contribution (sec Rcf. 16 for 
the parameters). We use a Fermi function for the ^°®Pb density as is given in 
Rcf. 24, while the same density is used for ^^O as that in Ref. 5. The double- 
folding potential roughly follows the Woods-Saxon potential with a=0.63 fm, 
especially at large distances, as would have been expected. It is in general 
much deeper than the Woods-Saxon potential, and accordingly some deviation 
from the Woods-Saxon potential is seen in fig. 2 at smaller distances. Figure 3 
compares fusion cross sections for the double-folding potential with those for 
the Woods-Saxon potentials. One can see that the double-folding potential 
leads to similar fusion cross sections to a Woods-Saxon potential with a=0.63 
fm. The double-folding potential alone, therefore, does not resolve the large 
diffuseness problem found in heavy-ion subbarrier fusion reactions. 
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Figure 4: Same as fig. 2, but for an error function potential. 



2.2 Complementary Error Function Potential 

We next consider the error function potential defined by 



VN{r) 



Vo 



erfc((r - Ro)/a), 



wliere tlie complementary error function is defined as 



erfc(x) = 1 j= 



dte 



(3) 



(4) 



It has been pointed out in Rcf. 11 that this potential falls more rapidly than 
the Woods-Saxon potential with increasing values of r while the curvature of 
the total potential remains similar. This is demonstrated in Fig. 4, where 
the error function potential with Vb=220 MeV, ro=0.92 fm, and a~2.9 fm 
is compared with the Woods-Saxon potential (note that the a parameter in 
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Figure 5: Same as fig. 3, but for an error function potential. 

the error function potential (3) does not have a unique correspondence to 
the surface diffuseness parameter in the Woods-Saxon potential) . At distances 
larger than r=10.5 fm, this potential is very close to the Woods-Saxon potential 
with a=0.63 fm, while at smaller distances it is similar to the Woods-Saxon 
potential with a=1.12 fm. As a consequence of the latter, the barrier curvature 
is nearly the same as that for the Woods-Saxon potential with a=1.12 fm, as 
can be seen in the upper panel of fig. 4. 

Fusion cross sections obtained with the error function potential are shown 
in fig. 5. It is interesting to notice that this potential gives fusion cross 
sections close to the Woods-Saxon potential with a=1.12 fm, thus providing 
an improved energy dependence of fusion cross sections. As is shown in fig. 
4, this potential is similar to the Woods-Saxon potential with a=0.63 fm at 
large distances, and therefore, would lead to similar scattering cross sections. 
The error function thus provides a promising potential which may describe the 
scattering and fusion processes simultaneously. 
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The main difference between the error function and the double folding 
potentials is that the former falls off as a function of inter-nuclear distance in 
a Gaussian way while the latter falls exponentially. In the folding procedure, 
it is implicitly assumed that the reaction takes place much faster than the 
dynamical density evolution. In this approximation, the projectile as well as 
the target densities remain unchanged even inside the Coulomb barrier and this 
leads to an exponential tail of the inter- nuclear potential E3 as a consequence 
of the fact that the nuclear density behaves exponentially at the surface. This 
approximation is referred to as the fmzen density approximation, and has been 
well tested in scattering processes. E3 However, if the error function potential 
reflected the "true "shape of the inter- nuclear potential, its Gaussian tail would 
suggest a deviation from the frozen density approximation (i.e., an exponential 
tail) for subbarrier fusion reactions, and effects of dynamical density evolution 
such as neck formationEZI would have to be taken into account. This is certainly 
an interesting problem, but is beyond our intention in this contribution and 
we leave it for future study. 



3 Effects of energy dissipation 

As an alternative origin for the large surface diffuseness problem, in this section 
we discuss the effects of energy dissipation. A very simple and phenomrun- 
logical way to mimic friction effects is to use the critical distance model. E3t3 
(See also Ref. 29). This model assumes that there is a certain distance Rc 
which has to be reached in order for fusion to occur. For each value of E, 
there is a corresponding critical angular momentum Ic above which Rc is in 
the classically forbidden region as angular momentum is increased. Only par- 
tial waves below Ic then contribute to fusion cross sections. Approximating the 
Coulomb barrier by a parabola, and limiting the partial wave summation up 
to the critical angular momentum, one_pbtains the expression for the fusion 
cross section in the Glas-Mosel model 113: 

<^fus{E) = —R, log I ^ _^ ^^2.iv,-E+iE-v.)R^JRl)/nn j ' (5) 

where fl, Rb and Vt are the curvature of the Coulomb barrier, the barrier 
position, and the barrier height, respectively. Vc is the value of the total 
potential at r = Rc- The critical angular momentum is given by 

Note that the Wong formula^ for the fusion cross section can be recovered by 
setting Rc to be infinity, that is equivalent to not restricting the partial wave 
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Figure 6: Fusion cross sections calculated with several models. The solid line is the predic- 
tion of the barrier penetration model which uses the Woods-Saxon potential with o=0.63 fm, 
while the dashed line is the result of the Glas-Mosel model. The parabolic approximation to 
the barrier penetration model is shown by the dotted line as a comparison. The thin solid 
line is the result of the surface friction model. 



summation in computing fusion cross sections. 

Figure 6 shows the results of analyses based on the Glas-Mosel model 
for the ^^O + ^"^Pb reaction. We first examine the validity of the parabolic 
approximation used in the Glas-Mosel model. The solid line is the result of 
the WS potential model with a=0.63 fm. If one approximates this poten- 
tial by a parabola, the barrier parameters are i?6=I1.94 fm, V{,=74.87 MeV, 
and ?if7=5.00 McV. The dotted line shows the fusion cross sections obtained 
with the Wong formula, which is based upon the parabolic approximation. 
Comparison between the solid and the dotted lines shows that the parabolic 
approximation works reasonably well in the energy range shown in the figure, 
although agreement is not perfect. 

We now vary i?c and independently to fit the experimental data, whilst 
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retaining the barrier parameters Rb, H, and hfl. The dashed hne is for 
i?c=10.39 fm and Vc—72A MeV. Ahhough the value of the critical distance Rc 
is larger than the systematics found in Ref. 19, this calculation reproduces the 
experimental data well at energies above the barrier, and the slope of fusion 
cross sections below the barrier is very similar to the prediction of the WS po- 
tential model with a=1.12 fm. Astonishingly, they are almost indistinguishable 
if they are plotted in the same figure at this scale. 

A more physical but still classical approach to incorporating dissipation 
is the surface friction model of Frobrich. Ej This model uses a single-folding 
potential, and friction coefficients which are propotional to the square of the 
first derivative of the nuclear potential. With slight adjustment of the inter- 
nuclear potential, and a changed tangential friction, quite good reproduction 
of the data was obtained (see thin solid line) . 

The surface-friction model is a macroscopic phenomenological model, and 
does not consider any particular origin of dissipation. One_pf the main origins 
would be inelastic excitations and/or particle transfer, til However, for the 
reaction ^^O + ^°®Pb, the effects of channel coupling are relatively weak due 
to the double closed shell structure of the projectile and the target. S Also it 
has been pointed out in Ref. 32 that energy dissipation due to particle transfer 
is largely compensated with the change of potential due to the polarization and 
thus the Coulomb barrier is reached almost without friction. Therefore neither 
collective excitations nor particle transfers would lead to the large dissipation 
effects, at least at the energies we are interested in. The only process which 
can be attributed to the origin of the dissipation would be a dynamical change 
in density of the reactants in the region where they overlap. This study again 
suggests a need to go beyond the frozen density approximation. 

4 Summary 

We studied possible origins of the recent finding that the experimental data 
for subbarrier fusion cross sections systematically favour a much larger surface 
diffuseness parameter for a Woods-Saxon form of inter-nuclear potential than 
the commonly accepted value of a ^0.63 fm. An obvious possibility would be 
that the inter-nuclear potential may deviate from the Woods-Saxon potential. 
We found that the double-folding potential roughly follows a Woods-Saxon 
potential with a=0.63 fm, and thus does not account for the large surface 
diffuseness parameter for subbarrier fusion. On the other hand, parameters 
could be adjusted for a complementary error function potential so that it looked 
like a WS with a large surface diffuseness at smaller distances while keeping 
a ~0.63 fm at the surface. The error function, therefore, provides a promising 
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potential which could simultaneously describe fusion and scattering. It falls 
off in a Gaussian way with increasing inter-nuclear distance whereas a WS 
potential falls off exponentially. The exponential fall-off is a consequence of the 
frozen density approximation. Hence, if the error function potential expresses 
the "true "shape of the inter- nuclear potential, this may suggest a dynamical 
evolution of density during fusion. As an alternative approach to the problem, 
we also examined the dissipation effects using the phenomenological critical 
distance model of Glas and Mosel as well as the surface friction model. The 
high energy data could be equally well reproduced both by these models and by 
using the WS potential with a large diffuseness parameter. The fact that the 
data could be explained both by changes in the potential, and the incorporation 
of friction, indicates that it may not be a simple matter to disentangle these two 
effects. We argue that the dynamical density change during the fusion process 
is one of the main origins of dissipation. This again suggests a deviation from 
the frozen density approximation, and thus a transition from the sudden to 
the adiabatic potentials. 123 In this connection, we would like to point out that 
in the adiabatic picture the inertia parameter may deviate from the reduced 
mass inside the barrier, as was discussed long ago by Mosel. cil Apparently, 
dynamical calculations are needed in order to draw a more conclusive picture 
for the surface diffuseness problem. Such calculations are now in progress, and 
we will report them in future publications. 
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